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Abstract. We analyze the Lyapunov exponents of U(l) gauge fields across the phase transition 
from the confinement to the Coulomb phase on the lattice which are initialized by quantum Monte 
Carlo simulations. We observe all features of a strange attractor with a tendency to regularity 
towards the continuum limit. Results are also displayed for the full spectrum of Lyapunov exponents 
of the SU(2) gauge system. 
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1. Classical chaotic dynamics from quantum Monte Carlo initial states. 

' Chaotic dynamics in general is characterized by the spectrum of Lyapunov exponents. 

C ■) , These exponents, if they are positive, reflect an exponential divergence of initially ad- 

QJ ' jacent configurations. In case of symmetries inherent in the Hamiltonian of the system 

there are corresponding zero values of these exponents. Finally negative exponents 
£f) ' belong to irrelevant directions in the phase space: perturbation components in these 

directions die out exponentially. Pure random gauge fields on the lattice show a char- 
acteristic Lyapunov spectrum consisting of one third of each kind of exponents [1] . 
The general definition of the Lyapunov exponent is based on a distance measure 
t^j- | d(t) in phase space, 

y—i • 

S '■ (1) L:= lim lim -In^r. 

g. t^cod(o)^oi d(0) 

, In case of conservative dynamics the sum of all Lyapunov exponents is zero according 

to Liouville's theorem, ^ L, ; = 0. For gauge field theories one utilizes the gauge in- 
variant distance measure consisting of the local differences of energy densities between 
two three-dimensional field configurations on the lattice: 
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(2) d:=^-J2 P \ trUp - tlU 'r\- 

Here the symbol J2 P stands for the sum over all Np plaquettes, so this distance 
is bound in the interval (0,27V) for the group SU(N). Up and U' P are the familiar 
plaquette variables, constructed from the basic link variables U Xi i, 

(3) U X)i = exp(aA Xti T) , 

located on lattice links pointing from the position x = {x\, X2,x?,) to x + aei. The gen- 
erator of the group U(l) is T = —ig and A x ^ is the vector potential. The elementary 
plaquette variable is constructed for a plaquette with a corner at x and lying in the 
ij-plane as U x ,ij = U Xt iU x+ ijU~l + j JJ X .. It is related to the magnetic field strength 
B x .k- 

(4) U X)i j = exp (e ljk a 2 B Xik T) . 
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The electric field strength E x ^ is related to the canonically conjugate momentum 
P x .i = aU x>i via 

(5) e x , = (vi, ,/:.,) . 

The Hamiltonian of the lattice gauge field system can be casted into the form 

(6) ff = E 

Here the scalar product stands for (A, B) = Hc{AB^). The staple variable V is a sum 
of triple products of elementary link variables closing a plaquette with the chosen link 
U. This way the Hamiltonian is formally written as a sum over link contributions and 
V plays the role of the classical force acting on the link variable U. 

We prepare the initial field configurations from a standard four-dimensional Eu- 
clidean Monte Carlo program on an N 3 x 4 lattice varying the inverse gauge coupling 
(3 oc g~ 2 [2]. We relate such four-dimensional Euclidean lattice field configurations to 
Minkowskian momenta and fields for the three-dimensional Hamiltonian simulation 
by selecting a fixed time slice of the four-dimensional lattice. 

2. Spectrum of the stability matrix. Instead of the classical determination 
of the Lyapunov exponent by the rescaling method outlined in the preceding section, 
we now use the monodromy matrix approach [3]. The Lyapunov spectrum Li is 
expressed in terms of the eigenvalues Aj of the monodromy matrix M: 

(7) L i= lim - f A^dt' , i = l,...,/ , 

t^oo t Jo 

where Aj(£) are the solutions of the characteristic equation for / degrees of freedom 

(8) det (Ai(t)t - M(t)) = 
at a given time t. Here M is the linear stability matrix, 

(9) M 

with the three-dimensional lattice action 53. 

Figure 1 displays the complex eigenvalues for selected U(l) gauge field configura- 
tions [4] prepared by a quantum Monte Carlo heat-bath algorithm. By investigating 
the time evolution of the monodromy matrix it turned out that the choice of the real 
time t does not affect the eigenvalues appreciably. In the confinement phase ((3 = 0.9) 
the eigenvalues lie on either the real or on the imaginary axes. This is a nice il- 
lustration of a strange attractor of a chaotic system. Positive Lyapunov exponents 
eject the trajectories from oscillating orbits provided by the imaginary eigenvalues. 
Negative Lyapunov exponents attract the trajectories keeping them confined in the 
basin. In the Coulomb phase (j3 = 1.1) the real Lyapunov exponents become rare to 
eventually vanish in the continuum limit [5]. In all cases the spectrum is symmetric 
with respect to the real and imaginary axes: the former property is due to the fact 
that the equations of motion are real, the latter is due to the Hamiltonian being con- 
servative (time independent). Also a number of zero- frequency modes occur which 



^(p,p) + i- 1 (u,v) 



( ad ati 
au ap 



\ au 



ap 
ap 



Full Spectrum of Lyapunov Exponents in Gauge Field Theory 



4 f 
3- 
2- 
1 ■ 

E ■ 

-1 ■ 
-2- 
-3- 

-4 1 



U(l) complex Lyapunov spectrum 



U(l) complex Lyapunov spectrum 





Re 



2 3 




Fig. 1. Full complex eigenvalue spectrum of the monodromy matrix for t = analyzing U(l) 
theory on a 4 3 lattice at /3 = 0.9 in the confinement (left) and at /3 = 1.1 in the Coulomb phase 
(right). 
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Fig. 2. Full complex eigenvalue spectrum of the monodromy matrix analyzing SU(2) theory on 
a 6 3 lattice at high energy g 2 aE = 0.8 (left) and at low energy g 2 aE = 0.1 (right) [3]. 



is connected to symmetry transformations commuting (in the Poisson bracket sense) 
with the Hamiltonian, such as time independent gauge transformations. 

Figure 2 depicts the complex eigenvalues for SU(2) gauge field configurations [3] 
taken from random starts with a high energy (g 2 aE — 0.8) and a low energy 
(g 2 aE = 0.1), respectively. Some caution is in order, since the gauge fields were 
not obtained from a Monte Carlo equilibrium configuration at definite coupling. Fur- 
ther, the entries of the monodromy matrix were constructed from derivatives with 
respect to components of the SU(2) matrix whereas the derivatives with respect to 
the phase were taken in the U(l) case. The SU(2) eigenvalues scatter in the complex 
plane but the symmetry of the spectrum with respect to the real and imaginary axes 
is still obvious. Again the prerequisites of a strange attractor are fulfilled and one 
observes some tendency to regularity at lower energy. 

Figure 3 shows the real part of the squared Lyapunov spectrum of U(l) gauge 
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Fig. 3. Real part of the squared spectrum of the stability matrix decreasing with increasing j3. 
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Fig. 4. Real part of the spectrum of the stability matrix with finite size scaling for increasing 



N [3]. 



fields for several couplings. Since the U(l) spectrum is purely real or purely imagi- 
nary, its square allows for a one-dimensional representation of its distribution. It is 
symmetric for (3 = and is shifted to negative values towards the Coulomb phase 
exhibiting an increasing number of zero modes with increasing ft. 

Figure 4 shows the real part of the Lyapunov spectrum of SU(2) gauge fields [3] 
extrapolated to 1/N — > from data taken at N = 2,3,4,5 and 6 at high energy 
(g 2 aE = 0.8). Since the SU(2) spectrum lies inside the complex plane, here only its 
real part is represented to gain a one-dimensional distribution. The overall pattern 
resembles that obtained earlier from smaller systems (N — 2, 3) with the rescaling 
method [6]. The structure of the ordered real part of the Lyapunov spectrum is 
similar at all energies considered, but the maximal point, aL max , scales with the 
energy, g 2 aE. 

3. Summary. This contribution concentrated on the spectrum of the mon- 
odromy matrix for (classical) compact U(l) theory. It exhibited an interplay between 
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positive, imaginary and negative Lyapunov exponents in the confinement phase chang- 
ing to a pure imaginary spectrum deep in the Coulomb phase leading to the regularity 
of the Maxwell theory. The spectrum for SU(2) theory was compared showing the 
characteristics of a strange attractor. Since the configurations were chosen from a 
random start corresponding to (3 = 0, it will be interesting to compute the spectrum 
of the stability matrix for finite (3 and to investigate if the non-Abelian theory stays 
chaotic in the continuum limit. 
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